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We compare the low-lying eigenmodes of the 0{a) improved Wilson-Dirac operator on quenched and dynamical 
configurations and investigate methods of probing the topological properties of gauge configurations. 


1. INTRODUCTION 


The eigenmodes of the Dirac operator carry in¬ 
formation about the topological content of the 
background gauge field. The Atiyah-Singer in¬ 
dex theorem tells us that in a sector with topo¬ 
logical charge Qt the Dirac operator has at least 
IQtl zero modes. In the presence of light dynam¬ 
ical fermions, the light quark determinant should 
suppress topological sectors with large \Qt\- As 
a consequence we expect a reduced value for the 
topological susceptibility 
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The topological charge is still difficult to deter¬ 
mine on the lattice. Various methods have been 
explored, usually gluonic methods. The latter, 
however, typically require smoothing procedures 
to reduce the ultraviolet fluctuations which do not 
take the fermionic part of the action, i.e. the pres¬ 
ence of dynamical fermions, into account. 

Fermionic methods, on the other hand, used to 
probe topological properties are expected to be 
sensitive to the effects caused by explicitly break¬ 
ing chiral symmetry when using Wilson fermions. 
These effects are potentially large on coarse lat- 
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tices, even when one reduces these effects using 
0(a) Symanzik improved fermions. 



Figure I. Eigenvalues of the massless improved 
Wilson-Dirac operator calculated on 0(150) 
quenched configurations (left) and dynamical 
(right) configurations. 


2. THE NON-HERMITIAN CASE 


The 0(a) improved Wilson-Dirac operator 


M{k) = 1 — k 
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is a non-hermitian operator. To calculate its 
eigenvalues \i we used the Arnold! algorithm [Q. 
To speed up convergence and to maximize the 
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Figure 2. The density of real modes (upper plot) 
and the topological susceptibility as a function of 
the largest real mode taken into account (lower 
plot) calculated on quenched configurations at 
/3 = 5.9. 



X 

Figure 3. Same as in Fig. || for quenched config¬ 
urations at /3 = 6 . 0 . 

number of real modes found by the algorithm, we 
transformed M using a “least-squares” polyno¬ 
mial of degree 40 1^ . This polynomial was tuned 
to make the eigenvalues with small real part lying 
in a band along the real axis fastest to converge. 
For our purpose it is crucial to increase the num¬ 
ber of calculated small real modes. We calculated 
0 ( 100 ) eigenvalues per configuration. 

The index theorem establishes a relation be¬ 
tween the difference in the number of real modes 
with negative and positive (n“) chirality 

Wi = (vi,^^Vi) and the topological charge: 

Qt = - n~. (3) 

However, on the lattice this relation can only 
hold approximately (if at all) in case of Wilson 
fermions. We can assume this approximation only 
to be reasonable if the spectrum of the Wilson- 
Dirac operator allows us to distinguish between 



Figure 4. Same as in Fig. for dynamical config¬ 
urations at (/3, K®®®-) = (5.2,0.1355). 

physical modes and unphysical modes due to the 
doublers. To check for this precondition we con¬ 
sidered the density of real modes. If our massless 
Dirac operator would be continuum-like we would 
expect the density of real modes being sharply 
peaked around zero and becoming small for larger 
modes. Comparing the densities at two differ¬ 
ent values of the gauge field coupling /3 in the 
quenched case we indeed see a more continuum¬ 
like behaviour for larger /3. In case of the dy¬ 
namical configurations we expect small eigenval¬ 
ues to be suppressed and therefore the density of 
real modes around zero to be smaller. This is in¬ 
deed what we find, see Figs. ||, 0 and ^ We now 
consider xt using Qt from Eq. (^ taking all real 
modes Ai < A into account. For larger values of 
the cut-off A one would expect Xt to become con¬ 
stant. In the quenched case Xt shows a reasonable 
plateau, which seems to improve for larger /?. The 
same analysis done on an ensemble of dynamical 
configurations does not show any plateau. 

3. THE HERMITIAN CASE 

The Wilson-Dirac operator becomes Hermitian 
when multiplied by 75 . To calculate 0(100) 
smallest eigenvalues of 75 M(k) we used again the 
Arnoldi algorithm with a Chebychev polynomial 
to speed up convergence of the eigenvalues with 
smallest modulus. Note that there is no trivial re¬ 
lation between the eigenvalues of 75 M(k) for dif¬ 
ferent values of k, while the eigenvalues of M{k) 
for different k are related by a trivial scale and 
shift operation. We used k = 0.1338 and 0.1353 
for /3 = 6.0 and (/3, k®®®-) = (5.2,0.1355), respec¬ 
tively, which corresponds to mps/my « 0.6 —0.7. 
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Table 1 

Simulation parameters. We used lattices of size 
V = 16^ X 32. The scale was set using tq = 0.5fm. 


p 

mps.sea [MeV] 

a [fm] 

#conf 

5.9 

quenched 

0.112 

0 ( 100 ) 

6.0 

quenched 

0.093 

0(150) 

5.2 

0.13550, 578(6) 

0.099 

0(150) 


The eigenvalues of the Hermitian operator can 
be used to calculate the topological charge Qt- 
Based on a chiral Ward identity one can define 

Qt = - lim TOqTr [ 75 M(TOq)“^] . (4) 

mq^O 

If the largest of the N calculated smallest eigen¬ 
values is large enough, we expect 

^ 1 

Tr[75M(TOq) 

to be a good approximation. For an extrapolation 
to iV = dim( 75 M) we performed a fit using the 
ansatz co -I- c_i |A|“^ ||]. 

While the theoretical basis for the definition 
of the topological charge using Eq. is much 
better, there are a number of disadvantages and 
open problems. For instance, calculations have 
to be performed at quark masses which are heavy 
enough to avoid effects from “exceptional config¬ 
urations” . Therefore, a chiral extrapolation is in 
principle required, but has not been done here. 
Initial results using smaller (and coarser) lattices 
show a very mild quark mass dependence |Q. Fi¬ 
nally, renormalization is a problem which has not 
yet been addressed. In Fig. H we compare the re¬ 
sults for xt calculated from the eigenvalues of M 
and 75 M. For quenched configurations at /3 = 6.0 
both methods give consistent results, indicating 
that the neglected renormalization factor is close 
to one. In the same figure we compare our results 
for Xt with those obtained in [|| using a gluonic 
method. The results agree surprisingly well. 

4. CONCLUSIONS 

We calculated the eigenvalues of the Wilson- 
Dirac operator M and the Hermitian operator 
75 M on dynamical and quenched gauge config¬ 
urations at similar lattice spacings. We found 
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Figure 5. Results for the topological susceptibil¬ 
ity using different fermionic methods (solid sym¬ 
bols, this work) and a gluonic method (open sym¬ 
bols, |]). 


clear signs for a surpression of small eigenmodes 
on configurations with dynamical fermions. We 
determined the topological susceptibility using 
two different definitions of the topological charge 
based on the index theorem and a chiral Ward 
identity. We found good agreement for quenched 
gauge configurations at /3 = 6.0 between the two 
methods. Our results indicate that the index 
theorem cannot hold on coarse lattices since the 
topological charge turns out to be ill-defined. Fi¬ 
nally, we found a surprisingly good agreement be¬ 
tween results for the topological susceptibility de¬ 
termined with gluonic and fermionic methods. 
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